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We study the dynamical behavior of disordered quantum-well-based semiconductor superlattices 
where the disorder is intentional and short-range correlated. We show that, whereas the transmission 
time of a particle grows exponentially with the number of wells in an usual disordered superlattice for 
any value of the incident particle energy, for specific values of the incident energy this time increases 
linearly when correlated disorder is included. As expected, those values of the energy coincide with 
a narrow subband of extended states predicted by the static calculations of Dominguez-Adame et 
al. [Phys. Rev. B 51, 14 359 (1994)]; such states are seen in our dynamical results to exhibit a 
ballistic regime, very close to the WKB approximation of a perfect superlattice. Fourier transform 
of the output signal for an incident Gaussian wave packet reveals a dramatic filtering of the original 
signal, which makes us confident that devices based on this property may be designed and used for 
nanotechnological applications. This is more so in view of the possibility of controlling the output 
band using a dc electric field, which we also discuss. In the conclusion we summarize our results 
and present an outlook for future developments arising from this work. 

PACS numbers; 73.20.Jc, 73.20.Dx, 71.20.-b, 85.42. -fm 



I. INTRODUCTION 

In recent years, there has been a growing interest in 
studies of disordered systems where the disorder presents 
some kind of correlation (see Ref. and references 
therein) . Aiming to find a physically realizable system of 
this type, Sanchez and Dominguez-Adame developed a 
simplified, continuous model in Ref. ^ for studying disor- 
dered semiconductor superlattices (SL's) where the disor- 
der exhibits short-range spatial correlations. In this par- 
ticular class of disordered SL's bands of extended states 
appear, opposite to the conventional view that in one- 
dimensional (ID) random systems almost all eigenstates 
are exponentially localized (see, e.g., Ref. |3|). Much more 
realistic calculations proved that these extended states 
are relevant to transport properties of actual superlat- 
tices, giving rise to large DC conductivities when the 
Fermi energy lies in one of these bands.u However, all 
those studies were carried out from a purely static view- 
point, and provided no information about the dynamics 
of electrons in this new type of nanostructures. 

In view of the lack of this kind of analysis, we un- 
dertook the study of the dynamical properties of elec- 
trons in these systems to complete the static picture, al- 
ready quite thorough. Thus, we compute the behavior 
of a wave packet incident on an intentionally disordered 
semiconductor SL by numerically solving the ID time- 
dependent Schrodinger equation for the complete Hamil- 
tonian (i.e., without analytical approximations) in the 
presence of an electric field. We explore several dynam- 



ical characteristics of our system, such as the tunneling 
times and the relation between the dwell time and the 
density of states.ETLl In addition, we estimate the charac- 
teristic time over which the resonant quasi-level can be 
established, showing that it is sufficiently large to allow 
the wave packet to tunnel close to the ballistic regime. 
We also consider the competition between quantum co- 
herence, preserved by correlated disorder, and the loss 
of quantum coherence due to an electric field acting on 
the SL. It^is important to clarify that loss of quantum 
coherencdj means in this context any elastic processes 
causing a complete localization of electronic states since 
we are not considering dissipative processes. Finally, we 
study the filter-like properties of these systems using the 
Fourier transform of the transmitted part of the wave 
packet and its dependence of the electric field, obtaining 
that it is possible to control the width and the center 
of the filtered band. It goes without saying that a cor- 
rect understanding of these properties is crucial from the 
perspective of technological applications of intentionally 
disordered SL's. 

The paper is organized as follows. In Sec. || v^KrPresent 



our model and summarize previous work of us,Qcl which 
we find convenient for a better understanding of the 
present paper, specifically as regards the behavior of the 
transmission coefficient, with and without electric field, 
for correlated and uncorrelated disordered SL's. The 



body of the paper is Sec. Ill where we present our dy- 
namical study of the system. We begin by examining the 
transmission probability and the transmission time for 
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the two different kinds of SL's. We compute the depen- 
dence of the transmission time with the size of the system 
in the WKB approximation for the balhstic regime and 
compare it with the numerical results. Most of the sec- 
tion is devoted to the relation between the mean dwell 
time and the density of states and, in addition, to the 
physical significance of the dwell time in this class of dis- 
ordered systems. We complete this characterization with 
a study of the spreading of the wave packet as a function 
of time. Following this equilibrium analysis, we devote 
Sec. 1^ to the study of the effects produced by the elec- 
tric field on the quantities presented in the last section, 
placing particular emphasis on the the filtering proper- 
ties of the correlated disordered SL's. Finally, in Sec. 0, 
we discuss our results and how can these be related to 
actual measurements to infer the main characteristics of 
the bands of the theoretically predicted extended states 
from experiments on SL's. We close the paper with a few 
prospects on future developments that may be attained 
starting from the present results. 



II. MODEL AND BACKGROUND 



have taken a = & = 32 A and a' = 26 A. The fraction of 
QWs of thickness a' is 40 % of the total number of QWs 
of the SL. This is not an essential parameter of the model 
as similar results are obtained taking other fractions. 

We now consider a single DQW as shown in Fig. |l| in 
an otherwise perfect and periodic SL. We showed ana- 
lytically in Ref. ^ that there is an specific energy value 
{Er) for which the so built SL is perfectly transparent, 
i.e. T{Er) — 1, where r is the transmission coefficient. 
The value of Er depends only on geometrical parameters 
(layer thicknesses) and it can be fixed at the fabrication 
stage. This result concerning resonant tunneling through 
a single DQW in an otherwise periodic SL does not im- 
ply that such a resonant phenomenon will survive in a 
disordered SL, that is, when more than one DQW's are 
randomly placed in th SL. The transfer-matrix formalism 
allows us to compute exactly, although not in a closed an- 
alytical fashion, the transmission coefficient in an arbi- 
trary SL. An example of the behavior of the transmission 
coefficient r around the resonant energy Er = 0.155 . . . 
is shown in Fig. ||(a) for a GaAs-Gao.65Alo.35As with 
N=200 barriers. 



We resume in this section previous results of usQtl for 
correlated disordered SL's in the stationary case, which 
will be useful for the discussion of the dynamical proper- 
ties which we address in the next section. For our present 
purposes, it is enough to focus on electron states close to 
the band gap with k|| = and use the one-band effective- 
mass framework to calculate the envelope-functions 



2m* dx^ 



+ VsLix) -eFx 



^Pix)=E^P{x), (1) 



where a explicit dependence of both E and ip{x) on quan- 
tum numbers is understood and they will be omitted in 
the rest of the paper. We have taken a constant effective- 
mass m* at the F valley although this is not a serious 
limitation as our description can be easily generalized to 
include two different effective masses. In the simplest pic- 
ture, the SL potential Vsl derives directly from the differ- 
ent energies of the conduction- and valence-band edges at 
the interfaces. A single quantum-well (QW) consists of a 
layer of thickness dj^ of a semiconductor A embedded in a 
semiconductor B. In our model of disordered SL, we con- 
sider that dA takes at random one of two values, a and a' . 
We call this a random SL (RSL). The thickness of layers 
B separating neighboring QWs is assumed to be the same 
in the whole SL, ds ~ b. A random dimer SL (DSL) is 
builiu by imposing the additional constraint that QWs 
of thickness a' appear only in pairs, called hereafter a 
dimer QW (DQW), as shown in Fig. |^. As a typical 
SL we have chosen a GaAs-Gao.65Alo.35As structure. In 
this case, the conduction-band offset is A£'c — 0.25 eV, 
and the effective mass is m* — 0.067 m, m being the 
electron mass. The origin of energies is taken at the 
GaAs conduction-band edge. In our computations we 



FIG. 1. Sdiematic diagram of tlie conduction-band profile 
of a SL containing a DQW. 

We next elucidated whether the physical mechanisms 
giving rise to delocalization in unperturbed systems are 
of relevance in the presence of an electric field, or the 
presence of the field destroyed the quantum coherence 
that exists at F = 0. To obtain the transmission co- 
efficient in the presence of an electric field, we develop 
a similar approach to that given in Ref. |l^. As usual in 
scattering problems, we assume an electron incident from 
the left and define the reflection, r, and transmission, i, 
amplitudes by the relationships 



ip{x) 



ikax 



re 



-ikax 



a; < 0, 
X > L, 



(2) 



where fcg 



2m*E/h, kl = {2m*{E + eFL)/h, and L 
is the length of the SL. The transmission coefficient is 
computed as t = {kL/ko)\t\'^ . Now we define ^p{x) — 
{\t\ ^/kL)q{x) ex.]3[i({x)], where q{x) and ({x) are real 
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functions. Inserting this factorization in Eq. (|l|) we have 
C.x{x) = q~'^{x) and 



2m* 



1 



VsLix)-eFx-E 



q{x) = 0. 
(3) 



This nonUnear differential equation must be supple- 
mented by appropriate boundary conditions. However, 
using Eq. (|^) this problem can be converted into a ini- 
tial conditions equation. In fact, it is straightforward to 
prove that 



(4) 



(5) 



q{L) = g,(L) = 0, 
and that the transmission coefficient is given by 
^ 4fcog^(0) 

^ l + 2fcog2(0) + fc2q4(0) + g2(0)<;2(0)- 

Hence, we can integrate numerically with initial con- 
ditions (^) backwards, from a; = L up to a; = 0, to obtain 
q{0) and ^^^(O), thus computing the transmission coeffi- 
cient for given incoming energy E and applied voltage 
V — FL. Figure ^(b) shows the transmission coefficient 
as a function of the incoming energy for a moderate value 
of the applied voltage F = lOkV/cm. We can see how 
the field shifts the mini-band to lower energies and de- 
stroy some of the quasi-bound states, but an important 
number of them survive. Then we have achieved the first 
goal of this paper: to show that the extended states that 
appear in DSL's survive in the presence of an electric 
field. Remembering that we proved previouslyu that this 
states also survive when interface roughness is taken into 
account, we can conclude that the delocalization due to 
structural correlations in the disorder is a very robust 
phenomena. In the next section we tackle the principal 
objective of this paper, namely to present a complete dy- 
namical study of the exciting properties of electrons in 
disordered DSL. 



(a) 



(b) 



Energy (eV) 

FIG. 2. Transmission coefficient r versus energy E 
for a DSL at (a) F = and (b) lOkV/cm. The 
GaAs-Gao.65Alo.35As SL consists of N = 200 barriers of 
6 = 32 A, whereas the thicknesses of QW are a = 32 A 
and a' = 26 A. 



III. DYNAMICAL RESULTS 



A. Numerical Method 



As we mentioned in the introduction, we are inter- 
ested in quantum diffusion of wave packets under an ap- 
plied electric field in semiconductor SL's. The equation 
which rules the evolution of the wave packet is the time- 
dependent Schrodinger equation. 



zn^^^^H{x)^{x,t), 



(6) 



where H{x) is the single-electron Hamiltonian given in 
(|l|). This equation has an elegant formal solution, given 
by 



(7) 



Using Cayley's form fcp,the finite difference representa- 
tion of the exponential 



l-jf-n{x)5t 



i + jf-n{x)5t 

we obtain the finite-difference equation 



1 



2fi 



nst ^ 



,fe+i 



(8) 



where we have replaced the wave function by its finite- 
difference approximation, in time (index k — 0,1,..., 
with — kSt) and in space (index j = 0, 1, . . . , A/") with 
Xj = j6x and Af the number of grid points). We will use 
a centered finite-difference approximation in x for 7i(x) 
and hence we have just a complex tridiagonal system. 
This method is commonly used in tbei solution of the 
time-dependent Schrodinger equationll3 because it en- 
sures strict norm conservation of the wave function at all 
times, and the error is only of the order {St'^). Norm con- 
servation has been used at every time step as a first test 
of the accuracy of results. We use a uniformly spaced set 
of spatial mesh points much larger than the SL's under 
consideration, and we transform the continuous bound- 
ary conditions, which read \E'(oo,t) — 'i>{—oo,t) = 0, to 
the corresponding discrete ones ^'^ij^ — — 0. Of 

course this approximation is valid only if we choose J\f 
sufficiently large to make sure than the wave function 
never comes close to the boundaries. We finally note 
that our initial wave function will be a superposition of 
plane waves of the form 



^'(a;,0) = [2Tr{Axf 



exp 



ikgx — {x — Xq)'^ 



A{Axy 



(9) 



where the average kinetic energy is E = h^kl/2m 



3 



To be published in Physical Review B 



^November 1996. 



B. Tunneling times and other dynamical tools 



The subject of tunneling, times is rich in contradictory 
definitions and results. BEl'tj When we measure the trans- 
mission time tT, we are trying to measure the time that 
a transmitted particle spent in the SL. The transmission 
time is straightforwardly obtained in the WKB limit for 
a ballistic electron, 




2 (A£;, - E) 



Xw {x) dx 



^ io V 2^ dx, (10) 

where Xb{z) and Xw{z) are the characteristic functions of 
the barriers and the wells, respectively. The mean dwell 
time tdw is 



tdw{E) = / dt \^{x,t)\'' dx, 



(11) 



and measures the average time spent by a wave packet 
in a given region of space. This time does not dis- 
tinguish between particles transmitted or reflected, and 
hence the mean dwell time becomes the transmission 
time of a transmitted particle when most of the wave 
packet is transmitted, as was pointed out by Biittiker 
and LandauerE 

Numerically, it is simple to measure tdw, and physi- 
cally is a powerful tool to measure the density of states, 
as can be shown thatD 



p{E) = \tdw{E). 
ttTl 



(12) 



According to Ref. ^ this relationship is only valid 
for symmetrical one-dimensional structures. For non- 
symmetrical structures it should be replaced by p{E) — 
2^ [^dwi^) ~^ ^dwi-^)] ' where the superscript refers to 
electrons coming from the right (r) or from the left (1). 
However we have found no differences between t'^^{E) 
and t^dwi-^) with the parameters we are using. 

Nevertheless, as Eq.(|l^) is only valid in a perfect bal- 
listic regime and the mean dwell time is only the trans- 
mission time in a idealized limit, we need to develop a 
method to measure tT- This method is based on the 
probability Pt that at time t the particle is found to 
have crossed the SL, 



Prit)^ 1^ mx,t)\'dx, 



(13) 



or the probability Pr that the particle is found to have 
been reflected back by the SL 



\ip{x, t)|^ dx, 



(14) 



To get an estimation of the spreading of the wave 
packet as a function of time we will use the time- 
dependent inverse participation ratio (IPR) and the 
mean-square displacement (cr), defined respectively as. 



IPR(t) = / \ij{x,t)fdx, (15a) 



with 



and will be explained in the next section. 



cr(i) = / {x^xy \iP{x,t)\'^ dx. (15b) 



x\ip{x,t)\^ dx. (16) 



Usually the IPR is a good estimation of the spatial ex- 
tent of electronic states. Delocalized states are expected 
to present small IPR (for long times IPR ~ ^/L), while 
localized states have larger IPR. The mean-square dis- 
placement is frequently also used to describe wave packet 
dynamics. In the asymptotic regime (t — > oo) one expects 
a behavior of the form (T{t) ~ P. Here < 7 < 1 for 
localization, 7=1 for ordinary diffusion, 1 < 7 < 2 for 
super-diffusion, and 7 = 2 for ballistic regime. The later 
is found in homogeneous systemaiJ. 



C. Quasi- ballistic scattering 

In this section we study the interaction of a Gaussian 
wave packet with average kinetic energy E, with the two 
different classes of disordered SL's, RSL and DSL, which 
we introduced in Sec. ||. For a RSL, we of course expect 
that the wave packet will be essentially reflected for any 
selected energy. However, in the case of a DSL we have 
two possible scenarios. On the one hand, if the dwell time 
is sufficiently large to allow a quasi-bound state of char- 
acteristic width r to be established, namely tdw — fi/T 
(see for example Ref. we expect that particles with 
energy close to the resonant one will be transmitted. If, 
on the contrary, the dwell time is not sufficiently large 
we never have a quasi-bound state and the behavior of 
the DSL will be the same that a RSL. A priori, we have 
no means to decide between these two possibilities, hence 
the necessity of the dynamical study we are summarizing 
here to clarify whether extended states do play a role in 
transport properties of DSL or not. 

Figure collects the results of a typical simulation of a 
wave packet for a DSL. In Fig. ||(a) we have a wave packet 
with a central energy oi E — 0.155 eV, very close to the 
resonant one obtained in Sec.||, traveling to impinge on 
a DSL. Some time afterwards, we can see in Fig. ^(b) 
that a small packet has emerged in the right part of the 
SL. We realize that the structure has filtered the initial 
wave packet, allowing only to pass the energies laying in 
the subband of extended states. We can confirm this in- 
terpretation by performing the Fourier transform of the 
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emergent wave packet and comparing it with the initial 
one as shown in Fig. |3|(c). We can see the emergent wave 
packet has an energy spectrum much narrower than the 
initial one, peaked around the resonance; this effect turns 
out to be much more dramatic the larger the SL is, but 
we preferred to keep within the limits of available super- 
lattices (note that = 50 in Fig. ||) instead of increasing 
the number of wells to get more spectacular results. 




(c) 



0.0 0.1 0.2 0.3 

Energy (eV) 

FIG. 3. The initial probability density 0)|'^, corre- 

sponding to a Gaussian function with Aa; = 200 A and energy 
E — Er — 0.155 eV, is shown in (a). The potential energy 
VsL(a;) is plotted as a function of x, for a DSL of A'' = 50 bar- 
riers of 6 = 32 A, whereas the thicknesses of QW are a = 32 
A and a' — 26 A. The inset in (a) is an enlarged view of a 
portion of the SL potential. The probability density 
at time t = 2 ps is shown in (b) . The Fourier transform for 
the initial wave packet (see (a)), \ip{E,0)\'^ , (dashed line) and 
for the transmitted packet through the DSL (see (b)) at time 
t = 2ps, \i)(E,t)\'^, (solid line) are shown in (c) as a function 
of the energy. 

We can understand better what is happening by look- 
ing at the dynamical evolution of the probability of trans- 
mission (reflection) PT{E,t) [Pi^{E,t)], i. e., the proba- 
bility of finding the particle in the right (left) side of the 
SL with energy E at time t. We notice that the station- 
ary transmission probability t{E) which we commented 
upon in Sec. |l| is just the limit of Pt{E, t), 



t{E) = C lim PT{E,t) 

t — >-oo 



(17) 



where C is a suitable normalization constant which de- 
pends on Ax and tends to unity as Aa; ^ cxi. In Fig- 
ure |[ we plot PT{E,t) and PR{E,t) as a function of 
time for the resonant energy E = 0.155 eV (solid line) 
and for E = 0.180 eV (dotted line) for a DSL, and for 
E = 0.155 eV (dashed line) for a RSL. For a RSL the 
results are similar for any energy; we have just selected 
Er as a typical behavior. For the DSL there is a great 
dependence of the energy. When we select an energy far 



from the resonant one we have a behavior similar to that 
the RSL. However, when we choose the resonant one, in a 
short interval of time Pt reaches practically its maximum 
value. In Fig. ||(b) we can see that this fast enhancement 
is due to the arrival of a compact packet corresponding 
to the components of the initial wave with energies closer 
to Er- Again, we have selected a small number of wells 
to allow an experimental verification of this results, and 
we have checked that the larger the number of wells is, 
the larger the differences between RSL and DSL. 



0.8 



0.6 




12 3 4 



dimers (out of resonance) 



Time (ps) 



FIG. 4. Dynamical transmission Pt (lower curves) and re- 
flection Pr (upper curves) probabilities are shown as a func- 
tion of time for a DSL with the same parameter as in Fig. 3, 
at energy E = Er = 0.155 eV (solid line) and 0.180 eV (dot- 
ted line) and for a RSL (with the same parameters that de 
DSL) at energy E = 0.155 eV (dashed line). The inset is an 
enlargement of the transmission probability for a RSL, the 
lowest curve in the plot. 

We are now going to use Pt and to find approx- 
imately the transmission time {tx)- We will choose the 
following convention: We obtain the time when the wave 
packet entries in the SL by finding the maximum value of 
dPii{t, E)/dt, which indicates the time when the proba- 
bility of finding the particle inside the SL is maximum, 
and we will fix our time origin tin at that instant (see 
Fig. H). This will be our time origin tin- As for the time 
when the wave packet exits the SL, we can think of it as 
the time when most of the particles transmitted are on 
the right part of the SL, i. e., the time tout when Pt is 
arriving to the final plateau (cf. Fig. ||). We obtain tout 
by finding the maximum value of dPxit, E)/dt- 

If the particle transmitted through the DSL is tunnel- 
ing through a ballistic channel induced by spatial cor- 
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relations in the disorder, the packet wiU pass the same 
amount of time in each well. Therefore, the time spent 
by the packet in passing through the whole SL would 
scale linearly with the number of wells, i. e., the length 
of the SL. One of the goals of this paper is to show not 
only that there is a significant enhancement of the trans- 
mission probability for particular values of the energy 
in DSL, but also that we are in the presence of a bal- 
listic transmission phenomenon in a disordered system, 
very close to the ideal WKB case for periodic SLs. This 
conclusion can be drawn from Fig. |. There we have 
plotted the transmission time as a function of the num- 
ber of wells, for both types of SL's, namely RSL and 
DSL, selecting an energy laying in the DSL mini-band. 
For comparison, we also show the results predicted by 
the WKB expression Eg. (|To|) . Remarkably enough, the 
DSL the behavior is purely linear and very close to that 
predicted by the WKB expression. On the contrary, for 
the RSL we have an exponential behavior characteris- 
tic of Anderson localized states. It thus becomes well 
established, from the dynamical view point, the nature 
of the DSL as a disordered system with good transport 
properties. 
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FIG. 5. We show the typical behavior of the (a) reflec- 
tion Pr and (b) transmission Pt probabilities versus time 
for the resonant energy in a DSL with the same parameters 
as in Fig. 3. The dashed line represents the derivatives (a) 
dPR{t,Er)/dt and (b) dPR{t, Er)/dt. The arrows mark the 
maximum values of the derivatives of Pr and Pt, which we 
take respectively as the initial time ti„ (i.e. when the wave 
packet goes inside the SL) and the outside time tout (i.e. when 
the the wave packet goes out of the SL). The transmission 
time is defined simply as tr = tout — tin- 



We now turn our attention to the mean dwell time. 
When we have a particle situated in an eigenstate, the 
mean dwell time in the DSL is exactly the transmission 
time Eujid has a clear relationship with the density of 
states.Q On the other hand, if we have a Gaussian wave 
packet the relation with the density of states is not at all 
evident. However, if we consider wave packets with small 
spread in energies we can expect that Eq. ^ continue to 
hold. In Fig. |^ we have plotted the tdw of an initial Gaus- 
sian wave packet of energy E = as a function of the 
number of wells. We can see that for the DSL the behav- 
ior is close to the linearity exhibited by the transmission 
time in Fig. |^ (solid line represents a linear fit). On the 
contrary, for the RSL the tdw exhibits a plateau because 
for this kind of SL's the dwell time is dominated by tn. 
In Fig. ^ we saw that for a RSL most part of the wave 
packet just penetrates in the SL a small number of wells; 
therefore tdw does not depend on the SL's size as soon as 
the SL is larger than those few wells. This result agrees 
with the typicai-consistency check for any definition of 
tunneling timejij where tdw is related to tx and tn by 
the expression. 



tdw = (1 - t) tji + TtT- 



(18) 



For a RSL when the number of wells grows, t goes to zero. 
In this case tdw is equal to ta, and hence the plateau ob- 
served in Fig. ^ 



& 4 



• RSL 

□ DSL 
WKB 



Number of wells 



FIG. 6. The transmission time tr at the resonant energy 
as a function of the number of wells for a DSL (squares) 
with a linear fit (dashed line) and a RSL (filled circles) with 
an exponential fit (solid line). Also depicted is the WKB ap- 
proximation for the DSL case (dot-dashed line), showing a 
good qualitative agreement (same linear behavior) with the 
numerical results for the DSL. Parameters for both SL's are 
the same as in Fig. 3. 

From the preceding considerations and results, we can 
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be confident that in the DSL, when the transmission 
and dwell times behave similarly and the propagation is 
quasi-ballistic, the aforementioned relationship between 
peaks in the density of states corresponding to peaks in 
the dwell time still holds, the only discrepancy being just 
a normalization constant, related to the amount of re- 
flected final states. Following this idea, we have plotted 
in Fig. H the density of states, obtained by using Eq. (12). 
We can see how we have a peak in tdw (and correspond- 
ingly in the density of states) for the resonant energy, 
as was expected for a packet that transmits through the 
whole SL by using the ballistic channel originated by the 
spatial correlation in a otherwise disordered SL. We note 
that this coincides with the predictions from the station- 
ary analysis in Ref. ^, thereby confirming again the re- 
sults in that paper. In particular, the shape of the density 
of states curve in Fig. |g is very similar to that obtained 
in our previous works. 



FIG. 7. Mean dwell time tdw at the resonant energy as a 
function of the number of wells A'^, for a DSL (circles) and 
a RSL (squares). A linear fit for the DSL (solid line) is also 
shown. Parameters for both SL's are the same as in Fig. 3. 

To conclude this zero field study of the dynamics of 
DSL's, we close the section by studying the spreading 
of the initial wave packet versus time for both kinds of 
SL's. In Fig. ^ we have plotted on a log-log scale the 
mean-square displacement cr as a function of time for a 
wave packet of energy E = Er incident in a DSL (solid 
line) and in a RSL (dot-dashed line), and for energy 
E — 0.180 eV impinging on a DSL (dashed line). At 
short times we can see a practically constant behavior of 
a which can be associated to the period while the wave 
packet is traveling towards the SL's. When the packet 
hits the SL, we see during a short time a decrease of cr, 
followed immediately by a rapid increasing of this mag- 
nitude. The decreasing is a consequence of the tails of 
the wave packet reaching the leading part, which is be- 
ing retained by its collision with the SL; once the whole 
packet is interacting with the SL, the behavior is close 
to power-like. We have fitted the results showed in Fig. 
^ for times larger than 1 ps to a function of the form 
a{t) ~ t"^ . We have obtained for the DSL at the reso- 
nant energy 7 = 1.120, i.e., we are in a super-diffusive 
regime, whereas for the DSL away from the resonant en- 
ergy 7 = 1.000, right at the limit between a localized 
regime and a ordinary diffusive one. Finally, for the RSL 
7 = 0.837, indicating that we are clearly in a localized 
regime. The a for the RSL is always much lower than for 
the DSL and increases with a larger approximate expo- 
nent, i. e., after some time (^ 1 ps) the packet is much 
more localized for the RSL than for the DSL (even out 
of the resonance). This is evidently a consequence of 
localization effects coming from the uncorrelated disor- 
der of the RSL, and whose influence is much less in the 
DSL case given the availability of extended states. Such 
a phenomenon further confirms the conclusions we have 



been drawing all along this section. We will come back 
to these results in the next section, when dealing with 
electric field effects. 



Energy (eV) 

FIG. 8. Density of states p as a function of the incident 
energy E obtained measuring tdw of a Gaussian wave packet 
of energy i5 in a DSL with the same parameters as in Fig. 3. 
Line is only a guide to the eye. 



Time (ps) 

FIG. 9. Mean-square displacement a of an initial Gaus- 
sian wave packet incident in a DSL at the resonance energy 
Er = 0.155 eV (solid line), and at energy E = 0.180 eV 
(dashed line) and for a RSL at energy E = 0.155 eV 
(dot-dashed line). Parameters for both SL's are the same 
as in Fig. 3. Note the log-log scale. 



IV. ELECTRIC FIELD EFFECTS 

The zero field simulations we have been presenting pro- 
vide an incomplete picture of electron dynamics in DSL's, 
as in this case it is obvious that technologically applica- 
ble phenomena would involve electric fields. Therefore, 
in this section we study the dependence of the dynamical 
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characteristics discussed in the previous paragraphs on 
an electric field. In Sec. || we explained that the sub- 
band of extended states appearing in DSL's is shifted 
to lower energies and reduces its width in the presence 
of moderate electric fields. We want to confirm that the 
correspondingly shifted quasi-bound states will have time 
enough to be established in the dynamical interaction of 
a Gaussian wave packet with a DSL potential when there 
is an applied field. To this end, in Fig. 10 we have plotted 
Pt for a Gaussian wave packet incident on a DSL with 
E — 0.155eV, for different values of the electric field F. 



sired values by selecting a particular value of the applied 
electric field. 




time (ps) 

FIG. 10. Transmission probability Pt as a function of time 
for a DSL with the same parameters as in Fig. 3, at energy 
E = Er = 0.155 eV, for different values of the electric field 
F = IkV/cm (dashed line) and 5kV/cm (dot-dashed line). 
For comparison we also show the result for F = Q (solid line). 



We can see that, at least for moderate fields, the dwell 
time is large enough {tdw — Tt-I^) to allow the quasi- 
bound states to be established thus permitting the trans- 
mission of the resonant components of the packet. The 
maximum value of Pt decreases with the field due to the 
shift of the mini-band and because the miniband becomes 
much narrower the larger the electric field applied is. In- 
terestingly, we can still find a dynamical resonant energy 
looking at the Fourier transform of the transmitted wave 
packet, as shown in Fig. |ll|. We want to stress that these 
results can be of interest for applications, because a DSL 
turns out to be a structure that works like an adaptive 
electronic filter, namely, by tuning properly the SL pa- 
rameters we can filter the energies contained in a narrow 
band. Moreover, this band can be displaced to the de- 
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FIG. 11. The Fourier transform for the transmitted packet 
versus energy in the presence of an electric field F = 10 kV/cm 
(dashed line); solid line shows the result at zero field. The in- 
sets show the transmission coefficient r as a function of energy 
for the two values of the electric field F for reference. The SL 
parameters are the same as in Fig. 3 



There is another aspect of the influences of electric 
fields on electron dynamics which is worth considering, 
namely the following: It is well known that, when an 
electric field is applied to a periodic SL, the localization 
of the initially extended states produces an oscillatory 
behavior of the wave packet, the so-called Bloch oscil- 
lations. Of course, Bloch oscillations require a quasi- 
perfect quantum coherence and a perfectly defined phase 
to be self-sustained in time. This is not the case in a DSL 
where electronic states increment their phase by a factor 
TT whenever they pass over a DQW.0 The question then 
arises as to what will the corresponding phenomenology 
in this case be. In order to answer this question, in Fig. 
|T^ we plot the IPR, defined in Eq. (15), as a function of 
time for (a) a perfect SL, (b) a DSL and (c) a RSL, in the 
presence of an electric field. The initial condition of these 
simulations was that, at t = 0, we placed a Gaussian wave 
packet with an energy oi E = 0.155 eV and Ax = 20 A in 
the center of each one of those SL's with 50 wells. In our 
case F = lOkV/cm, d = a -I- & = 64 A. In an ordered SL, 
the Bloch period will then be TbiocIi = h/eFd ^ 0.646 ps, 
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in perfect agreement with obtained from Fig. p^(a). For 
the disordered superlattices there is an oscillatory behav- 
ior at the beginning but in a short time the IPR achieves 
a randomly fluctuating, but stationary in mean, value. 
This indicates the existence of decoherence effects in both 
disordered lattices, the difference between the RSL and 
the DSL being that the latter shows a smaller mean value 
of the IPR, in agreement with the less localized character 
of its states. The remnants of oscillatory behavior for the 
disordered SL's are more clearly characterized by looking 
at the mean square displacement a as a function of the 
time, which is shown in Fig. |l3|. Again, we can see how in 
the DSL the wave packet is much more delocalized that 
in the RSL as a consequence of the presence of a narrow 
band of extended states. 
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FIG. 12. IPR of an initial Gaussian wave packet placed in 
a (a) perfect SL, (b) DSL, and (c) RSL as a function of time. 
For the three kind of SL's, the number of wells is A*' = 200 
and the rest of parameters are the same as in Fig. 3. 



FIG. 13. Mean-square displacement a of an initial Gaus- 
sian wave packet placed in a DSL (solid line) and RSL (dashed 
line), in the presence, of an electric field F — lOkV/cm, as a 
function of time. SL parameters are the same as in Fig. 12 



V. DISCUSSION AND CONCLUSIONS 

In this work we have successfully shown that the good 
transport properties predicted by previous static studies 
of SL's with correlated disorder give rise to correspond- 
ing dynamical phenomena of interest. To this end, we 
have reported on dynamical properties of electrons in 
intentionally disordered SL's, computed by using high- 
accuracy numerical methods to solve the time dependent 
Schrodinger equation for the complete Hamiltonian. In 
this respect, as the two main global conclusions of the 
present paper, we want to stress, first, that the dynami- 
cal results we show prove independently the existence of 
extended states with physical consequences in disordered 
systems, and second, that the validity of our previous 
static calculations in Refs. ^ and ^ to characterize elec- 
tron transmission through nanostructures has been set 
on firm grounds due to its perfect agreement with the 
dynamical analysis. 

Aside from the above general conclusion, which we 
draw from the consideration of a number of dynamical 
tools, we would like to summarize a few aspects more that 
we have learned from our simulation program. In partic- 
ular, we have proposed a method to find the transmission 
time by using the temporal transmission probability, pro- 
vided that such probability presents abrupt changes as a 
function of time. By means of this procedure to compute 
the transmission time, we have been able to show that 
the propagation of electrons with energies in the subband 
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of extended states of a DSL is ballistic, very similar to 
that of ordered SL's. This is a dramatic manifestation of 
delocalization by correlations, more so when compared to 
the exponential growth of the transmission time we have 
obtained for usual RSL's. In that regime, we have shown 
that the relationship between dwell times and density of 
states, holds for Gaussian wave packets, by computing 
the density of states and finding the same result as in 
our stationary calculations. Interestingly, the fact that 
correlations do not impede the phase decoherence of the 
wavepacket, the properties depending on symmetries of 
the system (translational invariance) are not recovered. 
This is the case, e. g., of Bloch oscillations. In any event, 
measurements of the IPR point out once more the differ- 
ences between DSL's and RSL's. All this characterization 
is confirmed by measuring the mean square displacement 
of electrons, which are seen to evolve faster in DSL's. 
Finally, we have also confirmed that low to moderate 
electric fields do not destroy the transport properties of 
DSL's, which is very important if DSL's are to be built 
and used for any practical purpose. 

To conclude, a few words are in order regarding possi- 
ble applications of the present work. It seems quite clear 
to us that several of our results can be useful for nan- 
otechnological devices with specific, special features. To 
begin with, the great difference of transmission times be- 
tween extended states and localized ones may provide a 
powerful tool for measure the extended character of the 
states in open systems. Besides, it can also be used to 
measure the amount and character of the disorder inher- 
ently present in any periodic SL, by obtaining the width 
of the band of extended states in the actual SL and com- 
paring it to the theoretically predicted one. However, 
what we think by far is the most promising application 
of DSL in nanotechnology has to do with their filter-like 
behavior. We have seen that it is possible, by means of an 
applied electric field, to control the center and width of 
the band of extended states, therefore allowing for a tun- 
able filtering of wavepackets, i. e., of electrons. This ca- 
pability, present already in practically achievable DSL's 
of some 50 wells, can be used to design a new family of 
electronic devices. In this respect, it is quite clear that a 
natural extension of this work would be to study the in- 
teraction of RSL and correlated disordered SL's, with an 
AC-electric field, using the complete Hamiltonian. Pre- 
liminary tight-binding resultalj appears to show exciting 
new phenomena in these structures. We envisage that 
appropriate choices of the frequency and/or intensity of 
the field can give rise to crucial changes in the filtering 



properties of DSL's. —Further work along these lines is 
currently in progress. E3 
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